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ABSTRACT 
It is found that Fe II emission contributes significantly to the optical and ultraviolet spectra of most active 
galactic nuclei. The origin of the optical/UV Fe II emission is still a question open to debate. The variability of 
Fe II would give clues to this origin. Using 7.5 yr spectroscopic monitoring data of one Palomer-Green (PG) 
quasi-stellar object (QSO), PG 1700+518, with strong optical Fe II emission, we obtain the light curves of the 
continuum /a (5 100) , Fe II , the broad component of H/3 , and the narrow component of H/3 by the spectral 
decomposition. Through the interpolation cross-correlation method, we calculate the time lags for light curves 
of Fe II , the total H/3 , the broad component of H/3 , and the narrow component of H/3 with respect to the 
continuum light curve. We find that the Fe II time lag in PG1700+518 is 209;!;J47 days, and the H/3 time lag 
cannot be determined. Assuming that Fe II and H/3 emission regions follow the virial relation between the 
time lag and the FWHM for the H/3 and Fe II emission lines, we can derive that the H/3 time lag is 148^[04 
days. The H/3 time lag calculated from the empirical luminosity-size relation is 222 days, which is consistent 
with our measured Fe II time lag. Considering the optical Fe II contribution, PG 1700+518 shares the same 
characteristic on the spectral slope variability as other 15 PG QSOs in our previous work, i.e., harder spectrum 
during brighter phase. 
Subject headings: black hole physics physics - galaxies: nuclei - quasars: emission lines 



1. INTRODUCTION 

The variability is a common phenomenon in quasi-stellar 
objects (QSOs) and provides a powerful constraint on their 
central engines. In the past two decades, the optical 
variability research focused on the spectral monitoring in- 
stead of the pure photometric monitoring. With the ac- 
tive galactic nuclei (AGNs) watch and the Palomer-Green 
(PG) QSOs spectrophotometrical monitoring projects, the 
reverberation mapping method, i.e., exploring the correla- 
tion between the emission lines and the continuum varia- 
tions, is used to investigate the inner st ructure in AGNs (e.g.. 



the Fe II emission between A4434 and A4684, strongly corre- 
lates with the so-called Eigenvector 1, which is su ggested to 
be driven by the accret ion rate (e.g., lBoroson & Gr een 1992J; 
Marziani et aTlbOOSal) . 

The origin of the optical/UV Fe II emission is still an open 
question. It is found that photoionized BLRs cannot produce 
the observed shape and strength of the optical Fe II emission 
and that the strength of UV Fe II cannot be explained unless 
considering the micro-turbulence of hundred s of km s~' or 
the co llisional exc itation in warm, dense g a s ([Baldwin et al. 



20041) . However, IVestergaard & Peterson I ( l2005b found the 



Blandford et al. 1982: Peterson 1993). It is found that mo 



tions of clo uds in the broad Une regions (BLRs) are viri - 
alized (e.g.. iKaspi et al. II2OO0I l2005UPeterson et al. Il2004 . 
With the line width of H/3 , Mg II , C IV from BLRs, the 
empirical size-luminosity relation derived from the mapping 
method is used to calculate the masses of their central su- 



permassive black holes (SMBHs; e.g.. iKaspi et al. I l2000t 



McLure & Jarvis 1 12004 iBian & Zhao Il2004 IPeterson et al' 



2004lGreene & Ho 11200.51)" 



It is found that the Fe II emission contributes significantly 
to the optical and ultraviolet spectra of most AGNs. Thou- 
sands of UV Fe II emission lines blend together to form a 
pseudocontinuum, resulting in the "small blue bump" around 
3000A wh en they are combin ed with Balmer continuum emis- 



correlation between the optical Fe II variance and the contin- 
uum variance and suggested that the optical Fe II is due to the 
line fluorescent in a photoionized plasma. It suggests that the 
optical Fe II line do not come frorn the s a me region as t he UV 
Fe II emission (e.g.. lKuehn et al. 1120081) . iMaoz et aU ( Il993h 
found that the reverberation time lag of UV Fe II in NGC 5548 
is about 10 days, similar to C IV time lag, smaller than the H/3 
time lag. The reverberation measurement for the optical Fe II 
emission has not fared so well. Some suggested that the opti- 
cal emission is produced in the same region as the other broad 
emission lines, and some suggested that it is in the outer por- 
tion of the BLRs because of narrower FWHM of Fe II with 



sion (e.g.. IWills et al. Ill985h . The optical Fe II would lead 



respect to H/3 (e.g..lLaor et al. II 19971: iMarziani et al. Il2003al : 



to two bumps in two s ides around the H/3 A4861A (e.g.. 



Boroson & Green 1992). It is found that the flux ratio of 



Fe II to H/3 , Rps, where the optical Fe II flux is the flux of 



Vestergaard & Peterson ll2005l:lKuehn et al. II2OO8I) . Recently. 
Hu et al. I (l2008allbl) did a systematic analysis of Fe II emis- 
sion in QSOs from flie Sloan Digital Sky Survey (SDSS), and 
found that the Fe 11 emission is redshifted with respect to the 
rest frame defined by the [O III] narrow emission line and H/3 
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Fig. 1. — Mean spectrum (the red curve) and rms spectrum (the green 
curve; multiphed by 8) for PG 1700+518. 

intermediate-width component is correlated with Fe II which 
locat es at the o uter portion of the BLRs. 

Ka spi et al. I POOO) gave the 7.5 yr spectroscopic moni- 
toring data for 17 PG QSOs. There is one PG QSO, PG 
1700- H518, with strongest optical Fe II emission an d Rpe = 
1.42 (iTurnshek et aTlligSSlilBoroson & Green1ll992h . Its H^ 



FWHM is 1846 ± 682km s"' jPeterson et al. Il2004 . and it is 



also called as a narrow line Seyfert 1 galaxy (NLSl). Using 
the Fe II template from one NLS 1 , 1 ZW 1 , we model the Fe II 
emission to investigate the Fe II variability and the relation to 
the continuum variability in PG 1 700H-5 1 8 . The data and anal- 
ysis are described in Section 2, the results are given in Section 
3, the discussion is given in Section 4, and the conclusions are 
presented in Section 5. All of the cosmological calculations in 
this paper assume Hq = 70km s~' Mpc"', i^M = 0.3, fl^ = 0.7. 

2. THE DATA AND ANALYSIS 

The spectroscopic monitoring data of PG 1700H-518 cover 
7.5 years from 1991 to 1998, which were done every 1-4 
months by using 2.3 m telescope at the Steward Observa- 
tory and 1 m telescope at the Wise Observatory. The total 
number of optical spectra for PG 1700H-518 is 39. The ob- 
servational wavelength coverage is from ^ 4000 to ~ 8000 A 
with a spectral resolution of ^ lOA. Spectra were calibrated 
to an absolute flux sc ale using siniultane ous observations of 



nearby standard stars (IKaspi et al. luOOOl) . The 39 spectra of 
PG 1700H-518 are available on the Web site '. 

In order to che ck its spectral variance, we calculate its m ean 
and rms spectra (iKaspi et al. Il2000t IPeterson et al. 1120041) . In 
Figure 1, the mean spectrum (the red curve) shows strong op- 
tical Fe II emission, and the rms spectrum (the green curve) 
shows variable emission for H/3 A4861, He I A5878, and Fe II 
(the blueward and redward of H/3 A4861). In the rms spec- 
trum, we can find Fe II features at ~ 4500A, 4924A and 
50 18 A, suggesting variable Fe II . The part of Fe II emis- 
sion between 4430A and 4770A in the rms spectrum is ob- 
vious than that between 5080 A and 5550A. It is due to the 
variable continuum slope, more variable in blueward of the 

' |http://wise-obs.tau.ac.il/~shai/PG/| 
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Fig. 2. — Example spectral decomposition for PG 1700+518. In the 
top panel, the black curve is the observed spectrum after the corrections of 
Galactic-extinction and the redshift, the red line is the sum of the power-law 
continuum and Fe II multiples (blue curves). The green ranges are our fit- 
ting windows. The bottom panel is the multi-Gaussian fits for the H/3 and [O 
III] lines. The red line is the sum of all multi-Gaussian (blue curves). The 
green curve is our fitting range of the pure H/3 and [O III] emissions after the 
subtraction of the power-law continuum and Fe II multiples. 



spectr um. It is consist ent with harder spectrum during bright 
phase (iPu et al. l2006l Figure 6). In the rms spectrum, we also 
find weak He II A4686. 

We use following steps to do the spectral decomposition, 
which have been us ed to analyze the spectra for a large QSOs 
sample from SDSS (IHu et al. IbOOSatlBian et al. 112008b . 

(1) First, the observed spectra are corrected for the Galactic 
extinction using Ay = 0. 116 from the NAS A/IPAC Extragalac- 
tic Database (NED), assuming an extinction curve of Cardelli, 
Clayton & IVlathis (1989; IR band) and O'Donnell (1994; op- 
tical band) with Ry = 3.1. Then the spectra are transformed 
into the rest frame by the redshift of 0.292. 

(2) The optical and ultraviolet Fe II template from the proto- 
type NLS 1 I ZW 1 is used to subtract the Fe II emission from 
the spectra ( IBoroson & Green 



pectra ( lljoroson & Ureen I ll992t IVestergaard & Wilkes 
200 lb . The I ZW 1 template is broadened by convolving with 
a Gaussian of various linewidths, the centroid wavelength 
shifts and fluxes. A power-law continuum is added in the 
fitting. The best modeling of the Fe II and the power-law 
continuum is found when x^ is minimized in the fitting win- 



Fe II variabiUty in PG QSO 1700+518 





22 


^ 


20 


Li- 
Li- 


18 
16 


„ 


16 


X 
li- 


14 
12 
10 




12 


„ 


10 



(i 



*** ft** iV 4 #**** **** 



:"^" I * \ii 



^^ \% 






tit* * 



**: *H 



Hp' 



^ 



ip 



^ J^* 



■I *' 



snv 



Mj *** **** 



*% 



t#* iH^ 



•*• 



A*^ 



Fig. 3.— Light curves of Fe II , H/3"+*, H/3* , H/3" , /a(5100) (from top to 
bottom). Emission-line fluxes are displayed in units of 10^'^erg s^'cm"^, the 
continuum flux determined from the power-law at 5100A is given in units of 
10""'ergs"'cm"2 A"'. 

dows: 4430-4770, 5080-5550A (see an example fit in Figure 
2). The monochromatic flux at 5100A, /\(5100) , is calculate 
from the power-law continuum. Because of the spectral cov- 
erage, we did not consider Balmer continuum (see the mean 
spectrum in Figure 1). 

(3) Considering weak [O III] AA4959,5007 lines, two sets 
of one Gaussian are used to model them. We take the same 
line width for each component, and fix the flux ratio of [O III] 
A4959 to [O III] A5007 to be 1:3. For the asymmetric profile 
of the H/3 profile, two-Gaussian is used to model the H/3 line, 
H/3^ and H/3" . The H/3^ and H/3" fluxes are calculated from 
integrating the corresponding fitting components. The flux for 
total H/3 , H/3"+\ is the sum of H/3* and H/3" fluxes. 

3. RESULT 

3.1. The light curves for the continuum, Fe II , ///3 

By IRAF-splot, the signal-to-noise ratios (S/Ns) between 
7400A and 7600A in the observational frame are measured 
for these 39 spectra. Two spectra (2nd and 35th spectra: 
pgl700910712.fits, pgl700980414.fits) are ig nored in our 



next analysis for their lower S/N less than 10 dKaspi et al 
2000l) . For the left 37 spectra, the distribution of S/N is 
25.3 ±9. The goodness of modeling of the Fe II and con- 
tinuum is tested by the elimination of Fe II features at A4924 
and A5017 (see Figure 2). We calculate the Fe II flux by in- 
tegrating the Fe II template fit between A4434 and A4684. In 
Figure 3, we show the fight curves of Fe II , H^"+^, H/3'' , H/3" 
, and /\(5100) (from top to bottom). 

We use the nor malized variability measure defined by 
Kaspi et al. I ( 120001) to compare the line variability to the con- 
tinuum variability, a^ = 100(cr^-(5^)'''^//, where / and a are 
the average and the rms of the flux in a given light curve, 
and 5 is the mean uncertainty in a given light curve. The a^ 
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Fig. 4. — Cross-correlation functions (CCFs) for the continuum-Fe 11 (top 
left), the continuum-H/3"^'' (top right), the continuum-H/i* (bottom left), and 
the continuum-H/3" (bottom right) for PG 1700-1-518. 
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Fig. 5. — Cross-correlation centroid distributions (CCCDs) for the 
continuum-Fe II cross-correlation (top left), the continuum-H/J cross- 
correlation (top right), the continuum-H/9 cross-correlation (bottom left), 
and the continuum-H/3" cross-coiTelation (bottom right) for PG 1700+518. 

are 6.73, 2.58, 8.46, 8.38, and 3.65 for the fight curves of 
/A(51 00),Fell,H/3\H /3" , and H/3"+*, respectively, in Fig- 
ure 3. lKaspi et al. I ( I2OOOI) showed that aj^ for the light curves 
of /\ (5 100) and H/3 are 6.8, 3.2, with respectively, which are 
consistent with our results. And the a^ of Fe II is smaller than 
that for /a (5 100) , H^* , H/3" , and H^"+^. 

3.2. Time lag from CCCD 

The measurement of the line time lag, t, is made by cross- 
correlating the emission line and c ontinuum light curve s. In 
order to compare with the result of dPeterson et al. 1120041) . we 
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use their code to do the line lag measurements. In Figure 
4, we give the interpolated cross-correlation functions (CCF) 
for the continuum-Fe II (top left), the continuum-H/?""*"* (top 
right), the continuum-H/3''(bottom left), and the continuum- 
H/3" (bottom right) for PG 1700+518. We find a peak in the 
CCF for Fe II light curve, which can be used to determine the 
time lag for Fe II curve (Figures 4, 5). Because of many peaks 
with positive lag times and/or peaks with negative lag times 
in the CCFs for other three curves, we can't determine their 
time lags (Figures 4 and 5). 

Through the Monte Carlo FR/RSS method, the uncertainty 
of T can be determined (Peterson et al. 2004). In the code, we 
adopt a minimum correlation coefficient of 0.4, the centroid 
threshold of 0.8 ^max, and the number of trials is 3000 (Peter- 
son et al. 2004). In Figure 5, we give the cross-correlation 
centroid distributions (CCCDs) for the continuum-Fe II 
cross-correlation (top left) and the continuum-H/?""^* cross- 
correlation (top right) for PG 1700+518, as well as that for 
H/3'' and H^" (bottom panels). It is obvious that the CCCD for 
Fe II hav e a narrow positive pe ak. Following the suggestion 
given by ( iPeterson et al. 1120041) . we calculate the mean of the 
CCCD in all valid trials as the Fe II centroid time lag, as well 
as its upper and lower uncertainties. We find that Fe II time 
lag ,TFeii, in PG 1700+518 is 270!j|[J days. Its mean CCF /?n,ax 
is 0.54 ±0.08. There are many positive peaks and/or strong 
negative peaks for H/3 , H/3" , and H/S* . Therefore, we cannot 
give the fine lags for H/3 , H/?" , and H/S* . With the redshift of 
0.292 for PG 1700+8 15, in the rest frame, the Fe II time lag in 
PG1700+518 is 209![*;^' days. Because the H^ time lag cannot 
be determined, we do not know whether the region emitting 
Fe II is located outside of the regio n emitting t he broad H/3 
lines (Marzianiet al. 2003bt Vester gaard & Pe terson 2005 : 
Popovic et al. I|20071IHu et al. Il2008bl: iKuehn et al. Il2008h . 



4. DISCUSSION 



4.1. 



The Fe U fitting method 

In the analysis of the optical Fe II light curve for PG 
1700+518, we fit the optical spectrum by the Fe II template 
instead of directly calculating the optical Fe II flu x in the se- 



lected wavelength range (e.g.. lKuehn et al. 1120081) . In model 



ing Fe II emission, we simultaneo usly model the pow er-law 



continuum, which is different from lWang et al. I ( 120051) in the 



Fe II analysis of an NLS 1 NGC 405 1 . We also consider vari- 
ous FWHMs, centroid wavelength shifts and fluxes in the con- 
volving of the Fe II template. The accuracy of the measure- 
ment for the continuum shape depends on the wavelength cov- 
erage. For these 37 spectra of PG 1700+518, the wavelength 
covera ge is mainly between ^ 3500Aand ^ 6000A in the rest 



frame. IVestergaard & Peterson I (120051) found that the optical 
Fe II feature to the blue of H/3 is contaminated by He I A4471 
and He II A4686 fines. For PG 1700+518, the He I A4471 and 
He II A4686 lines are not strong (Figure 1). Therefore, we use 
file fitting windows of 4430-4770 and 5080-5550A to exclude 
the emission fines of H^ , H7 A4340, and [O III] AA4959, 
5007 (Figure 2). When we mask the region of He II A4686 in 
the fitting, the fitting result is almost the same. We tried Fe II 



template of lVeron-Cettv et al. I (120041) . it is almost the same to 
thatoflZWl. 

4.2. The optical Fe II emitting region 



Kuehn et al. I ( 120081) presented the reverberation analysis of 
optical Fe II for Ark 120. They gave the light curves of the 
blue/red side of Fe II , H/3 , and continuum by setting the 
measurement windows (Figure 1 in Kuehn et al. 2008). Al- 
though the Fe II cross-correlation function is very broad and 
flat-topped, they suggested that the optical Fe II -emitting re- 
gion, ^ 320 days, is several tim es larger than the H/3 zone 



{^ 57 davs) lKuehn et al. I (120081) found that it is difficult to 
constrain the FWHM of optical Fe II for Ark 120 because 
of its very smooth Fe II emission (Figure 8 in Kuehn et al. 
2008). Modeling the Fe II emission in PG 1700+518, we find 
that the mean value of Fe II FWHM is 1554 ± 1 10 km s"' . 
Because the change of H/3 profile due to the Fe II contribu- 
tion is not too much (Figure 1), we adopted the H/3 FW HM 
value of 1846 ± 682 km s"' by IPeterson et aTl (|2004 . We 
find that (FWHMHa/FWHMpeii)^ = 1.41. Assuming that Fe 
II and H/3 emission regions follow the virial relation between 
the time lag and the FWHM for the H/3 and Fe II emission 
lines, we can derive that the H/3 time lag is 148;!;ig4 days. We 
also find taht the new estimate d H/3 time lag i s cons istent with 

their 



(|2009[ 



^blr-ALa(5100) relation bv iBentz et al. 
Figure 5). 

Considerin g the host contribution in /a(5100) , 
Bentz et al. I ( 120091) suggested a new relation be- 
tween BLRs size and ALa(5100) , log/?BLR = 
(-21.3^?) + (0.51 9!!!°g^)logALx(5100) (It 

the 
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days). 
(I2OOOI) gave the average flux of /a (5 100) 
(between 6520 and 6570 A) without excluding Fe II contri- 
bution, (22.0 ± 1.5) X 10"'*erg s'^cm'^A"'. After excluding 
Fe II contribution, we find that the value of /a (5 100) at 
5100A is (21.4 ± 1.5) X lO-'^'ergs-'cm-^A"'. Therefore, 
the Fe 11 correction is very small for /a(5100 ) . Cor rected 
the contribution from starlight, iBentz et al. I (120091) gave 
/a(5100) as (18.5 ± 1.5)-"'erg s-icm-^A"' and ALa(5100) as 
3.63 X lO^^erg s"' (the starlight contribution is about 16% in 
its to tal flux). The expec ted 7?blr from /?blr- ALa(5 100) rela- 
tion dBentz et al. Il2009t) and ALa(5100) of 3.63 x lO^^erg s"' 
is 222 It-days. Considering the larger uncertainty of intercept 
in this relation (about 3 in log/?BLR), this result is consistent 
with our estimated H/3 time lag, 148 
2OO9I see their Figure 5). 



^tuM days ( IBentz et al. 



If we take the FWHM/time lag uncertainties into consid- 
eration, the Fe II emission region is located near the H/3 
emission region, not conclusively located outside of the H/3 
emission region. Kuehn et al. (2008) suggested that op- 
tical Fe II emission is possibly produced at the dust subli- 
mation radius, T^dust = 476 x [Lboi/lO'^^erg s"']°-^ (It-days) 



dEfitzur & Shlosmanl I2OO6I) . By Lboi = 9ALa(5100), for 



PG 1700+5 18, we find that the dust sublimation radius /?dust ^ 
2868 (It -days), which is much larger than the radius indi- 
cated by the Fe II emission lag time. 



Kaspi et al. I (12000) measured the H/3 flux between 6 120 A 



and 6410 A in the observational frame (also in Peterson 
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Fig. 6. — Spectral slope, a, vs. /x(5100) , fy oc v^°' . Considering the 
error of a, the red line is the best linear fit. The green Hne is the best linear 
fit excluding two discrete red points. 
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tion. 
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2004). Their H /3 fluxes include Fe II contribu- 



Kaspi et al. I (I2OOOI) gave H/3 flux as (18.88 ±0.99) x 
erg s~'cm"^. Removing Fe 11 contamination in this re- 
gion, the H/3 flux is (14.22 ± 1.01) x lO'^^g^g ^-icvcT-. With 
our new light curve for H/3 , we cannot determine the time lag 
for H/3 emission line, and we can determine the time lag for 
Fe II line. It is possible due to: (l)H/3 line has a asymmetric 
profile (Figure 2; narrow and broad components), suggesting 
that H/3 is coming from the region with a very broad size and 
Fe II is coming from the region with a narrow size . There- 



fore, w e ca n detect Fe II time lag . The H/3 flux i riKaspi et al 



teOOOh andlPeterson et al. I ( |2004 includes Fe II contribution. 
(2 iPeterson et al. I (120041) found two peaks in its CCCD (see 
their Figure 15) and suggested that the peak at zero is due to 
correlated error. The spectral S/Ns are not high and the time 
sampling in the light curves is not good. More data and higher 
S/N spectra are needed in the future. 
With respect to the H/3 time lag of 252 days suggested by 



Peterson et al. I (120041) . smaller estimated H/3 time lag of 148 
days, which leads to the smaller black hole mass estimation in 
the logarithm, is decreased by 0.23 dex. However, consider- 
ing the uncertainties of time lag and the mass calculation, our 
results are consistent with that from lPeterson et al. I (|2004|) . 



4.3. The relation between the spectral index and /a(5 100) 

In our previous work lPu et al. I (120061) . we investigated the 
relation between the spectral index and /x(5100) and found 
almost all (15/17) PG QSOs showed an anti-correlation be- 
tween them, except PG 1700H-518 and PG 1229H-204. For PG 
1229H-204, it is due to the system difference in two telescopes. 
For PG 1700H-518, we suggested that it is maybe due to the Fe 
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II contribution dPu et al. l2006h . Here, we give this relation for 
PG 1700H-518 when Fe II contribution is carefully removed. 
We find that there is a strong anti-correlation between them 
(see Figure 6). The spearman coefficient R is -0.33, with a 
probability of pauw < 0.05 for rejecting the null hypothesis 
of no correlation. If two discrete red points are excluded, R 
is 0.55, and p^uw < 5.8 x 10"^. Therefore, after considering 
Fe 11 contribution, PG 1700H-518 shares the same characteris- 
tic on spectral s lope variability as other 15 PG QSOs in our 
previous paper (IPu et al. 1120061) . i.e., harder spectrum during 
brighter phase (Hubeny et al. 2000). 

5. CONCLUSION 

With the spectral decomposition of 39 spectra of PG 
1700H-518 with the strong Fe II emission, we investigate the 
Fe 11 variability and its time lag. The main conclusions can be 
summarized as follows: (1) we give light curves of /x(5100) 
, Fe II , H/3'' , H/3" , and H/3 "+*, as well as the mean and 
rms spectra for PG 1700H-518. With the normalized variabil- 
ity measure, a^, we find that all components are variable. (2) 
With the code of Peterson et al. (2004), we find that Fe II time 



lag in PG1700H-518 is 209![57 days, and H/3 time lag cannot 
be determined. (3) Considering the uncertainties of time lags, 
the expected H/3 time lag from the empirical luminosity-size 
relation is 221.6 It-days, consistent with our measured Fe 11 
time lag. If we take FWHM/time lag uncertainties into con- 
sideration, Fe II emission region is located near the H/3 emis- 
sion region, not conclusively located outside of the H/3 emis- 
sion region. (4) Assuming that Fe 11 and H/3 emission regions 
follow the virial relation between the time lag and the FWHM 
for the H/3 and Fe II emission lines, we can derive that the 
H/3 time lag is 148;!;[Q4 days. With respect to the H/3 time lag 
of 252 days suggested by Peterson et al. (2004), smaller H/3 
time lag, which leads to the black hole mass estimation in the 
logarithm, is decreased by 0.23 dex. (5) After considering 
Fe 11 contribution, PG 1700H-518 shares the same characteris- 
tic on spectral slope variabihty to other 15 PG QSOs in our 
previous work (IPu et al. 1120061) . i.e., harder spectrum during 
brighter phase. 
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